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We present an experimental study of Z-scan measurements of the nonlinear response of cold Cs 
atoms in a magneto-optical trap. Numerical simulations of the Z-scan signal accounting for the 
nonuniform atomic density agree with experimental results at large probe beam detunings. At 
small probe detunings a spatially varying radiation force modifies the Z-scan signal. We also show 
that the measured scan is sensitive to the atomic polarization distribution. 
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I. INTRODUCTION 

Clouds of cold atomic vapors produced by laser cool- 
ing and trapping techniques provide strongly nonlinear 
effects when probed with near resonant light beams. 
Soon after the first demonstration of a magneto optical 
trap(MOT) 1] it was recognized that the cold and dense 
clouds were useful for studies of nonlinear spectroscopy 

generation of non-classical light 4] and other nonlinear 
effects 5, 6j. In this work we study nonlinear focusing of 
a near resonant probe beam taking into account the spa- 
tial dependence of the atomic density and the effect of the 
probe induced radiative forces on the MOT distribution. 

From the perspective of nonlinear optics a cold vapor 
in a MOT is an interesting alternative to the nonlinear- 
ity of a hot vapor in a cell. A typical MOT provides a 
peak density of order 10^ — 10^^ cm~^, and an interac- 
tion length of a few mm. On the other hand a heated 
vapor cell can easily have atomic densities of 10^^ cm~^ 
and interaction lengths of 10 cm, or more, so that much 
stronger nonlinearities can be achieved in a vapor cell 
than in a MOT. Nonetheless there are a number of rea- 
sons for looking more closely at the optical nonlinearity 
provided by a cold vapor. To begin with the relative 
strength of the MOT nonlinearity compared with a hot 
cell is larger than the above estimates would indicate be- 
cause of Doppler broadening effects. In the presence of 
the large inhomogeneous Doppler width in a hot cell it is 
typical to detune the probe beam by a GHz or more to 
avoid strong absorption. Conversely the Doppler limited 
linewidth of cold atoms in a MOT is of order the homoge- 
neous linewidth or less so that much smaller detunings of 
order tens of MHz can be used. The full saturated non- 
linearity can therefore be achieved with a much weaker 
probe beam in a MOT than in a vapor cell. Addition- 
ally, the large detunings used in vapor cells imply that 
detailed modeling of nonlinear effects in Alkali vapors 
must account for the plethora of hyperfine and Zeeman 
levels. Such modeling, as has been done by several groups 
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0, -) involves numerical integration of hundreds of cou- 
pled equations for the density matrix elements. Only in 
special cases, using e.g. a buffer gas to create large pres- 
sure broadening, can a simplified two-level type model 
provide an accurate description of nonlinear propagation 
effects 9]. Working in a MOT with a small probe beam 
detuning the probe interacts strongly with only a single 
atomic transition so that the main effects of nonlinear 
propagation can be achieved using a compact two-level 
atomic model. As we show below this is only partially 
true since radiative forces and population distribution 
among Zeeman sublevels lead to observable features in 
the measured probe transmission. 

In addition a cold atomic vapor potentially provides 
a qualitatively different nonlinearity than a hot vapor 
does because the mechanical effects of light can result 
in strong modification of the atomic density distribu- 
tion, which in turn feeds back on the nonlinearity seen 
by the optical beam. Indeed some experiments already 
observed reshaping of the MOT cloud due to radiation 
trapping effects 0, E3 While such density redis- 

tribution effects in both position and momentum space 
may also occur in hot vapors, as in the collective atomic 
recoil lase r ITp. IT^l H^ . and radiation pressure induced 
dispersion T5J, such effects are potentially much more 
pronounced in cold vapors where momentum transfer 
from the light beams is significant. In particular we ex- 
pect that complex spatial structures can be formed due 
to coupled instabilities of the light and density distri- 
butions in a fashion analogous to the effects that have 
been predicted for light interacting with coherent mat- 
ter waves 1^ Some recent related work has shown 
evidence of nonlinear focusing in a MOT 18], and possi- 
bly structure formation in experiments that include the 
effects of cavity feedback [Toj. In Sec. |^we discuss the 
relevance of the present measurements in the context of 
observation of coupled light and matter instabilities. 

Our primary interest in the present work is a detailed 
study of nonlinear focusing and defocusing of a tightly fo- 
cused probe beam that propagates through a MOT. The 
Z-scan technique was originally developed 20j for charac- 
terization of thin samples of nonlinear materials. Since 
the MOT cloud is localized to a region of a few mm in 



2 



thickness we can easily apply this technique for character- 
ization of the MOT nonlinearity. The theoretical frame- 
work based on a two-level model is described in Sec. ^ 
Z-scan measurements were taken with a Cs MOT using 
the procedures discussed in Sec. IIIII The experimental 
and theoretical results are compared in Sec. I^where we 
also compare additional measurements and calculations 
of reshaping of the transverse profile of the probe beam 
after propagation through the MOT. 

II. NUMERICAL MODEL 



(a) Top view 




In Z-scan measurements, the transverse profile of a 
laser beam passing through a nonlinear sample is investi- 
gated. In the presence of self-focusing or self-defocusing 
the transmittance through a small aperture placed after 
the medium exhibits a S-shaped dependence on the po- 
sition of the beam waist with respect to the nonlinear 
sample. Follo wing the original demonstration of the Z- 
scan technique^, 21] there have been a large number of 
theoretical studies of the Z-scan method that take into 
account different types of nonlinear response and con- 
sider different characteristic ratios between the Rayleigh 
length of the probe beam and the thickness of the non- 
linear sample |2ll (2^, [2^ [2^, 12^, 123 . The interaction of 
a probe beam with a cloud of cold two-level atoms is de- 
scribed by a saturable mixed absorptive and dispersive 
nonlinearity with a susceptibility of the form j2 71] 



X(r) na{r)Wo 
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Here na{r) is the density of atoms at position r, Wq is 
the population difference in thermal equilibrium, 7 is the 
homogeneous linewidth, A = a; — 0;^ is the difference be- 
tween the probe frequency co = 27rc/A and the atomic 
transition frequency uJa^ A is the wavelength of the laser 
beam in vacuum, Ig is the on resonance saturation inten- 
sity, and / is the optical intensity. 

None of the existing theoretical treatments can be di- 
rectly used for our situation. Although the work of 
Bian, et al.|2^ studied the Z-scan behavior for a sat- 
urable nonlinearity, it was restricted to no absorption 
and weak saturation. In work to be published elsewhere 
we derive analytical expressions for the Z-scan curve for 
a susceptibility of the type given in Eq. How- 
ever, we find that in order to obtain good agreement 
with experimental measurements it is necessary to take 
account of the spatial variation of the density in the 
MOT cloud. We have therefore relied on direct numer- 
ical simulations to compare with experimental results. 
To do so we assume a scalar probe beam of the form 
E = [A{T)/2]e<^^-^^^ + [A*(r)/2]e-<^^-^^) and invoke 
the paraxial and slowly varying envelope approximations 
to arrive at the wave equation 
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(b) Timing sequence 
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FIG. 1: (color on-line) Z-scan experimental setup (a) and 
timing sequence (b). D and D^e/ are signal and reference de- 
tectors respectively and D fi is the fluorescence detector used 
to monitor the number of cold atoms. 



where V^^ = jdx^ + jdy^ with x^y the transverse 
coordinates and k = 27r/A. 

Equation Q was solved numerically on a 128 x 128 or 
256 X 256 point transverse grid using a split-step spectral 
code. Propagation from the atomic cloud to the pinhole 
plane was calculated by solving Q with the right hand 
side equal to zero. As the experimental free propaga- 
tion distance of 15 cm resulted in the beam spilling over 
the edges of the computational window, calculations were 
done with a distance of 3 cm, and a pinhole diameter of 
1/5 the actual size. The numerical solutions obtained 
in this way describe the interaction of the probe beam 
with an unperturbed MOT. Since there is no integration 
over the Doppler profile of the atoms we are implicitly 
assuming they are at rest. In reality the atoms are accel- 
erated due to the absorption of photons from the beam, 
so the laser frequency seen by the atoms is shifted to 
the red. A model including this Doppler effect should be 
used when the probe beam interacts with the cloud for a 
time corresponding to many absorption/emission cycles. 
We discuss the implications of the radiative force for the 
Z-scan curves in Sec. IIVI 



III. EXPERIMENTAL SETUP 

The experimental setup is shown in Fig. ^a). A stan- 
dard MOT was loaded directly from a background Cs va- 
por. The trapping and repumping beams were obtained 
from two external cavity diode laser systems. The trap- 
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ping beams were detuned from the F = 4 ^ = 5 cy- 
cling transition by A ~ —27, where 7/(27r) = 5.2 MHz. 
The Gaussian diameter of each trapping beam was 2.5 
cm with a peak intensity of 2.4 mW/cm^, and the beams 
were retroreflected after passage through the MOT. The 
peak saturation parameter for all six beams was 13.1. 
The repumping beam was tuned to the F = 3 ^ F^ = 4 
transition. The magnetic field gradient was 9 G/cm along 
the vertical y direction while the probe beam propagated 
horizontally along z. Time of flight measurements gave a 
typical MOT temperature of 60 /iK. 

In order to model the Z-scan data accurately some care 
was taken to characterize in detail the spatial distribu- 
tion of trapped atoms. Fluorescence measurements of 
the MOT cloud taken with a camera placed on the x 
axis revealed a flattened density profile indicative of ra- 
diation trapping effects |3, 0, [2^. This type of profile 
has previously been modeled with a Fermi-Dirac type 
distribution 29]. We chose to use an expansion with more 
fitting parameters of the form 

_2zi _2z± _2z£ 

, , a^e ^bze ^ c^e 
na{y,z) = UaO ^— 

+ ^2 + Cz 

_ 2y^ _ 2^4 _ 2^6 
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^ Trrr^ — 

CLy \ Oy ~r Cy 

where n^o is the peak density, Wy are the Gaussian 
radii along z and y and the a, 6, c are fit parameters along 
the two axes. The trapped Cs atoms formed an ellipsoidal 
cloud, which was modeled as a density distribution Eq. 
with typical Gaussian diameters of 2wz = 7.3 mm 
and 2wy = 6.0 mm as shown in Fig. [21 The fit residuals, 
as seen in the lower plots in Fig. [3 were small in the 
center of the MOT and reached 20% at the very edges of 
the cloud. Since the axis of the B field coils was along y 
we assumed that the x and z density profiles were equal. 
As the size of the probe beam was small compared to the 
width of the cloud we approximated na(r) by na{z) given 
by the first line of Eq. in our numerical simulations. 

The total number of trapped atoms was 2.5 x 10^ as 
measured by an optical pumping method ^ . This 
number was also measured and monitored from day to 
day by recording the fluorescence signal using a lens and 
calibrated photodetector 28]. The result from the flu- 
orescence signal, taking into account the corrections for 
the atomic polarizability distribution discussed in Ref. 
[28] was about 2.5 times lower than that from the opti- 
cal pumping measurement. The number of cold atoms 
measured by the fluorescence method varied by up to 
15% from day to day. The peak atomic density deter- 
mined from the number measurement using the optical 
pumping method, and the size of the cloud using the 
data shown in Fig. 13 was typically 5.5 x 10^^ cm~^, and 
2.2 X 10^^ cm~^ using the fluorescence measurement of 
the number of atoms. We believe the optical pumping 
method to be more accurate since it does not rely on any 
assumptions about the polarization of the MOT cloud. 
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FIG. 2: (color on line) Intensity profiles (dots) of the Cs 
atomic cloud along z (a) and y (b) directions, with curve 
fitting results (solid lines). Frames (c) and (d) show the cor- 
responding fitting errors. 



The Z-scan probe beam was derived from the trapping 
laser and frequency shifted to the desired detuning with 
an acousto-optic modulator. The beam was then spa- 
tially filtered with an optical fiber before focusing into 
the MOT with a lens fscan = 400 mm mounted on a 
translation stage. The Gaussian radius of the beam at 
the focus was wq = 24.5 /im. A 1.0 mm diameter pin- 
hole was placed 15 cm away from the center of the MOT 
(outside the vacuum chamber). The transmitted field 
through the pinhole was measured by photodetector D. 
To measure the transmittance of the probe beam, the 
trapping beams and the probe beam were turned on se- 
quentially, as shown in Fig. ^b). During the measure- 
ment, the trapping beams were turned off first and the 
repumping beam was left on to pump the atoms into 
the F = 4 state. The transmittance of the pinhole at 
different times after the probe beam was turned on was 
measured as D{t)/Drei' 

In the experiment, the lens was scanned instead of the 
nonlinear sample as shown in Fig. ^a). Compared with 
traditional Z-scan measurements, there are two points to 
consider. First, note that a movement of the lens in the 
-\-z direction is equivalent to a movement of the cloud in 
the —z direction. Therefore, Z-scan curves obtained in 
this experiment have the opposite configuration to the 
traditional ones, e.g. a peak followed by a valley shows 
a self-focusing nonlinearity. Second, since the lens was 
scanned, the position of the beam waist changed, which 
affected the linear (without cold atoms) transmittance. 
To take this effect into account, we recorded 
with the MOT on and off (by turning the magnetic field 
on and off) so that the normalized Z-scan curve was given 
by as a function of z. 

Before discussing the Z-scan measurements we note 
that transmission scans with the pinhole removed can be 
used for calibration of the peak atomic density. Figure 
01 shows the measured transmittance curve with the pin- 
hole removed and the probe beam tuned on resonance 
with the transition F = 4 ^ = 5. The peak inten- 
sity was 1.07 W/cm^ giving a peak saturation parameter 
of 975. Measurements were taken for both circular (a^ 
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and a~) and linear polarized beams. The measurements 
show that the probe beam was strongly absorbed when 
the beam waist was far from the center of the atomic 
cloud. As the beam waist was moved closer to the cloud 
center the intensity of the beam increased and the ab- 
sorption saturated, so that the transmittance increased. 
The results using right and left hand circular polarization 
were very similar to each other, while the curve for the 
linear beam is broader than those with circular beams 
and its peak transmittance is a little lower. The solid 
line is the numerical calculation under the experimental 
conditions using a two- level atom model with I5 = 1.1 
mW/cm^, and the dashed line is the calculation with Is 
= 1.6 mW/cm^. Thus the solid line corresponds to a cir- 
cularly polarized probe with the assumption of complete 
optical pumping of the atoms into the uif = ±4 levels 
(I5 = 1.1 mW/cm^), while the dashed line corresponds to 
linear polarization with the assumption of uniform distri- 
bution among the Zeeman sublevels (Ig = 1.6 mW/cm^). 

The best fit to the data implies a peak atomic density 
of Uao = 2.3 X 10^^ cm~^. This is lower than the measure- 
ment based on the optical pumping method, and slightly 
higher than the measurement using the fluorescence sig- 
nal. As can be seen from the figure, the numerical result 
for a circularly polarized beam agrees reasonably with the 
data, although the width of the calculated scan is about 
15% narrower than the data. The agreement between 
calculation and experiment for linear probe polarization 
is slightly worse. 



IV. EXPERIMENTAL RESULTS 

A. Z-Scan measurements 

Figures El and [SI show Z-scans measured for self- 
defocusing and self-focusing nonlinearities. Measure- 




z (mm) 

FIG. 3: (color on-line) Measured on resonance transmission 
curve of a Cs cloud at A = 852 nm 6/is after the probe beam 
was turned on with different polarizations. The solid line is 
the calculated result for a circularly polarized beam and the 
dashed line is for a linearly polarized beam. The Gaussian 
diameter of the atomic cloud was 2wz =6.2 mm and — 
0.85,5^ = 0.03, = 0.07. 



ments were taken with right hand circular (RHC) and left 
hand circular (LHC) probe beam polarizations at each 
detuning. In all the experiments, the peak intensity of 
the probe beam was 1.07 W/cm^. The solid lines show 
the numerically calculated curves which assumed com- 
plete optical pumping of the atoms into the tuf = ±4 lev- 
els as in the transmission measurements described above. 

Concentrating first on the self-defocusing case with 
RHC polarization shown in Fig. ^ a,c we see that at 
short times after turning off the MOT the shape of the 
curves is similar to that calculated numerically, while the 
peak density which was used as a fitting parameter was 
about half as large as that measured using the methods 
described in Sec. IIIII We believe that the discrepancy 
is due to partial atomic polarization as discussed below. 
The curves have a pronounced peak when the probe waist 
is positioned past the center of the cloud at z ~ 7 mm and 
only a small dip, which almost disappears as the detun- 
ing is decreased from —47 to —37. The departure from 
the traditional S shaped form that is obtained for a pure 
Kerr nonlinearity 20j , is due to the nonlinear absorption. 
The measurements were taken in a regime of strong sat- 
uration at the probe beam waist, sa = 26.4 at A = —37 
and sa = 15.0 at A = —47 where the saturation param- 
eter is SA(r) = I{t)/Isa and Is a = /s(l + 4A^/7^) is the 
saturation intensity at finite detuning ^27^. As we move 
to larger detuning the saturation parameter decreases, 
so the effects of nonlinear absorption are diminished and 
the Z-scan curves tend towards the normal symmetric S 
shape. 

At long times of t = 2.6 ms the MOT cloud has 
partially dispersed and the Z-scan curve is flatter and 
broader with an almost complete disappearance of the 
valley at negative z. In the intermediate regime of t = 
112 /is the peak at z = 7 mm increases by 5-10% com- 
pared to the value at t = 16 /is. This increase is consistent 
with focusing, and a local increase in density, due to the 
radiation force 

7 '^(^) 
Fsp(r) = ?i k - ^ ^ _ ^ v(r)) 2/72 + s(r) ' ^"^^ 

where ^(r) = I(y)/Is and v(r) is the atomic velocity. 
As shown in Ref. 32] a radiation force that decreases 
with distance into the MOT due to absorption of the 
pushing beam results in focusing of the atoms after a 
finite time. Near the focal plane of the probe beam where 
the saturation parameter is large the peak light induced 
acceleration is amax = fikj/{2m)^ with m the atomic 
m ass. The t ime to be pushed a distance 5z is thus 5t — 
\/25z/am3.^' which for ^^^Cs evaluates to 5t :^ 130 /is 
for 5z = 0.5 mm. This is consistent with the observation 
of a larger Z-scan signal at a probe interaction time of 
t 112 /is. 

Turning now to the self-focusing case with LHC polar- 
ization shown in Fig. [nia,c we see similar results as in the 
self-defocusing case at t = 16 /is except that the trans- 
mittance peak now appears for z negative as expected. 
The valley in the transmittance is also more pronounced 
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FIG. 4: Measured Z-scan of a Cs cloud at A = 852 nm at different time delays with detunings A = —4j (a,b) and A = —87 (c,d) 
with respect to the F = 4 ^ = 5 transition. The data marked by open circles, triangles and crosses show the transmittance 
at t = 16 /iS, 112 fis and 2.6 ms after the probe beam was turned on respectively. The MOT parameters used in the calculation 
were 2wz = 7.2 mm, az = 0.32, bz = 0.27, Cz = 0.26 (measured from camera images), and riao = 1.0 x 10^° cm~^ (fitting 
parameter) . 
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FIG. 5: Measured Z-scan of a Cs cloud at different time delays with detunings A = +47 (a,b) and A = +87 (c,d) with 
respect to the F = 4 ^ F' = 5 transition. The data marked by open circles, triangles and crosses show the transmittance at 
t = 16 //s, 112 //s and 2.6 ms after the probe beam was turned on respectively. The MOT parameters used in the calculation 
2wz = 7.4 mm, az = —0.88, bz = 2.2, Cz = —1.1 (measured from camera images), and Uao = 0.9 x 10^° cm~^ (fitting parameter). 



than for self-defocusing particularly at the larger detun- 
ing of A = +47. At the intermediate time of t = 112 /is 
we see an enhanced peak which again is consistent with 
focusing due to the radiation pressure force. 



However at the longest time of t = 2.6 ms we see the 
overall flattening of the Z-scan curve together with an 
unexpected secondary peak in the transmission that ap- 
pears for z positive. This secondary peak can be quali- 



6 



tatively explained using Eq. (Q. The radiation pressure 
force accelerates the atoms and Doppler shifts them to 
the red. Although the initial detuning is positive for self- 
focusing, after some time a fraction of the atoms will 
be Doppler shifted to an effective negative detuning and 
will act to defocus the light which gives a transmission 
peak at positive z. We can estimate the time for this to 
occur by noting that the peak time dependent Doppler 
shift is kv(i) = ka^^yi. A Doppler shift of 7 is reached at 
= 2m /{hk'^) which evaluates to = 78 jas for -"^^^Cs. 
Thus for A = 37 and t > we expect to see some 
evidence of self-defocusing, which is consistent with the 
experimental data in Fig. [51 a, c at t = 2.6 ms. The effect 
is smaller when we go to A = +47 since the scatter- 
ing force is weaker and the Doppler shift must be larger 
in order to change the sign of the detuning. Note that 
no secondary features appear in Fig. 21 where we start 
with red detuning since the Doppler shifts only move the 
atoms even further out of resonance. 

The measurements discussed above were repeated with 
opposite helicity of the probe beam as shown in the 
right hand columns of Figs. ^ and [5J The strength of 
the Z-scan signal was substantially different for the two 
probe beam helicities. For example for red detuning at 
A = —47 the peak to valley Z-scan signal at t = 16 /is 
from Fig. 21 a, b was 0.57 for RHC polarization but only 
0.39 for LHC polarization. For blue detuning at A = 47 
the peak to valley Z-scan signal at t = 16 /as from Fig. 
El a,b was 0.57 for LHC polarization but only 0.37 for 
RHC polarization. Thus for red detuning the strongest 
effect was obtained with RHC polarization which was 
the opposite polarization as the copropagating (j+ trap- 
ping beams shown in Fig. d while for blue detuning 
the strongest effect was obtained with LHC polarization 
which was the same polarization as the copropagating 
(7+ trapping beams shown in Fig. ^ We considered sev- 
eral possible reasons for this unexpected dependence on 
probe beam helicity. All measurements were taken with 
the MOT magnetic field still on. However, the Zeeman 
shift of the cycling transition at the edge of the MOT 
cloud in the x — z plane was less than 3 MHz which is 
not large enough to explain any polarization or Zeeman 
dependence. 

The numerical simulations were done with a two-level 
model which implicitly assumes complete pumping of the 
atomic population to the lower level of the cycling transi- 
tion. Before application of the probe beam the atomic po- 
larization will be distributed across Zeeman levels. The 
spatial localization provided by a MOT is due to the 
fact that (j+ and a~ polarizations interact more strongly 
with the atoms on different sides of the cloud. This nat- 
urally leads to a spatial variation of the atomic polariza- 
tion. Indeed a rate equation model 33] predicts a polar- 
ization of 20% or more at the edges of a Rb MOT. At 
very short times after the probe beam is turned on, or 
if the optical pumping due to the probe beam were only 
partially successful, we would expect a polarization de- 
pendent interaction that depended on the helicity of the 



probe. However, Figs. 21 and [51 show that the helicity 
of the probe that interacts most strongly with the MOT 
switches when we change the sign of the detuning. 

This effect can be explained qualitatively in the fol- 
lowing way. Looking at Fig. Elwe see that the effective 
off-resonance saturation parameter of the probe beam 
varies by many orders of magnitude across the atomic 
cloud. The volume that determines the nonlinear diffrac- 
tion of the probe beam can reasonably be taken to extend 
out to off-axis distances of a beam waist where the off- 
resonance saturation parameter is as small as 0.001 at 
the edges of the cloud for a centered probe beam. Let 
us assume the atomic population before application of 
the probe beam is distributed across the Zeeman lev- 
els. Then the longest optical pumping time, that for 
transferring an atom in mi? = —4 to mi? = 4, will be, 
neglecting the excited state branching ratios, tpump ^ 

(2/7)[(l + 5A)/sA]E™::i4l/|Cn^..™.+lP where the 
Clebsch-Gordan coefficients are normalized such that 

CmF=4.,mp,^5 = L We find tpump ^ 20 /iS at 5A = 1 

and tpump ~ 1 ms at 5a = 0.01. The implication of this 
estimate is that the probe beam does not completely po- 
larize the volume of the MOT that it interacts with, even 
at times as long as several hundred jas. At the latest time 
shown in Figs. 21 and O t = 2.6 ms, the optical pumping 
is substantially complete, and there is only a very small 
difference between the data taken with opposite probe 
beam helicities, except for the effects of the radiation 
force as discussed above. 

The Z-scan signal is intrinsically due to the nonlinear 
effects of self-focusing and self-defocusing. These effects 
are strongest when a circularly polarized probe interacts 
with a fully polarized atomic sample. However, the radi- 
ation forces that shift the detuning to the red, are also 
maximized for a fully polarized atomic sample. When 
the probe is blue detuned the radiation forces will ini- 
tially increase the strength of the interaction so that a 
probe polarization that couples strongly to the atomic 
polarization will give a larger effect. In the opposite case 
of red detuning of the probe the pushing forces only serve 
to decrease the strength of the interaction so that the Z- 
scan signal will be strongest when the pushing is reduced. 
Although the MOT has opposite atomic polarization for 
positive and negative z and is not expected to have any 
net polarization imbalance averaged over the cloud, the 
probe beam is attenuated as it propagates so that the 
front edge of the cloud has a stronger impact on the Z- 
scan signal. We thus expect a larger signal when a red de- 
tuned probe has a polarization that gives a weaker atomic 
coupling, and less pushing forces, on the front side of the 
cloud, and a larger signal when a blue detuned probe 
has a polarization that gives a stronger atomic coupling 
on the front side of the cloud. These arguments suggest 
that a red-detuned probe will give a larger signal when 
it has the opposite helicity of the trapping beam that 
has the same momentum projection along i, and that a 
blue-detuned probe will give a larger signal when it has 
the same helicity as the trapping beam. This is indeed 
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what is observed in Figs. 2|and[5J 

As a further check on this explanation we changed the 
sign of the magnetic field and reversed the helicities of all 
the trapping laser beams. We then found that the data 
were the same as that measured previously provided we 
also fiipped the helicity of the probe beam. This sup- 
ports the conclusion that the dependence on probe beam 
polarization is due to the spatial distribution of atomic 
polarization inside the MOT in the presence of imperfect 
optical pumping and radiative forces. It is interesting 
to note that the number density deduced from transmis- 
sion scans shown in Fig. |31 would be higher, and thus 
in closer agreement with the other density measurements 
discussed in Sec. IIIII if we took into account a partial po- 
larization of the atoms. However the transmission scans, 
which were taken without a pinhole, and therefore de- 
pend only on the total transmission, and not the shape 
of the wave front, show no helicity dependence. We con- 
clude that Z-scan measurements provide a signal that is 
comparatively sensitive to the atomic polarization distri- 
bution. The extent to which a Z-scan could be used to 
measure the polarization distribution quantitatively re- 
mains an open question. 




-8 -6 -4 -2 2 4 6 8 
z (mm) 

FIG. 6: Saturation parameter of the probe beam in the atomic 
cloud at different radial positions p — yx^-T^. w{z) — 
wq^I + X^z'^ / [ti'^Wq) is the z-dependent probe beam waist. 



B. Spatial redistribution of intensity 

In addition to the Z-scan transmittance measurements, 
the transmitted probe beam far field distribution was 
observed directly using a CCD camera. To do so, the 
pinhole and photodetector D in Fig. Q were removed, 
and a lens (/ = 100 mm) was used to image a plane 
at z = 72 mm after the cloud center onto the camera. 
Pictures were recorded with the beam waist at differ- 
ent positions relative to the center of the cloud for self- 
defocusing and self-focusing cases as shown in Figs. [7| 
andlHl To illustrate the self-defocusing or self-focusing 
of the beam, the difference between the intensity distri- 
butions with and without cold atoms is shown. In Fig. 
[7|we see the result with the beam waist near the center 
of the cloud {z = —1 mm) and at the edge of the cloud 



close to the detector {z = +6 mm). At z = — 1 mm, the 
beam gets focused faster due to the nonlinearity, so the 
far field transmitted beam has an additional divergence 
compared to the linear case. Thus the central part of the 
picture is dark. At z = +6 mm, the beam becomes less 
focused due to the same effect, so that in the far field the 
beam is more converged relative to the linear case, which 
gives a bright region in the center of the picture. For the 
two pictures in Fig. [SI we have similar results with the 
roles of positive and negative z interchanged due to the 
opposite sign of the nonlinearity. 

The lower parts of Figs. [7| and |H1 show that numeri- 
cal calculations predict the transverse beam profiles with 
an accuracy similar to that seen for the Z-scan trans- 
mittance curves. We see that the agreement is best for a 
self-defocusing nonlinearity while in the self-focusing case 
there is a tendency towards localized maxima and min- 
ima in the transverse profiles. Similar phenomena were 
reported in Ref. [l^. Additional numerical calculations 
at 5-10 times higher densities reveal strong filament ation 
of the transmitted beam. Future work will investigate 
propagation effects in this regime experimentally. 



V. DISCUSSION 

We have used Z-scan measurements to characterize 
the nonlinear optical response of a cold atomic cloud. 
We show that the measured data agree with a two-level 
atomic model at short times provided the probe helicity 
is chosen to match the helicity of the trapping beams on 
the front side of the cloud. At longer times modifications 
to the Z-scan occur because of radiation pressure forces. 
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FIG. 7: Transmitted probe beam intensity for A = —37- 
The pictures show the intensity minus the intensity without 
a MOT for z — —1 mm and +6 mm. The bottom figures 
show the line profiles across the center of the beam, which are 
normalized to the center intensity of the beam without cold 
atoms. The dashed line is the result of numerical calculations. 
The peak MOT density was riao = 0.7 x 10^° cm"^ while 
Wz,az,bz,Cz were the same as in Fig. |5| 
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X (mm) X (mm) 

FIG. 8: Transmitted probe beam intensity for A = +87 at 
z — —3 mm and +2 mm. All other parameters are the same 
as in Fig. 



This results in some additional focusing of the cloud and 
the appearance of a secondary Z-scan peak when a blue 
detuned probe accelerates the atoms past resonance, to 
give a partially red detuned response. Transverse profiles 
of the probe beam show focusing and defocusing features 
consistent with the Z-scan transmittance measurements. 
Future work will study clouds with larger optical thick- 
ness where we expect modulational instabilities of the 
probe beam to lead to small scale modifications of the 
atomic distribution. 
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